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Abstract 

The off mass shell continuation of dual amplitude with Mandelstam analyticity (DAMA) is pro- 
posed. The modified DAMA (M-DAMA) preserves all the attractive properties of DAMA, such as 
its pole structure and Regge asymptotics, and leads to a generalized dual amplitude A{s,t,Q^). In 
such a way we complete a unified "two-dimensionally dual" picture of strong interaction m^lOllll- 
This generalized amplitude can be checked in the known kinematical limits, i.e. it should reduce to 
' the ordinary dual amplitude on mass shell, and to the nuclear structure function when t — 0. We 

' fix the Q^-dependence in M-DAMA by comparing the structure function F2, resulting from it, with 

phenomenological parameterizations. The results of M-DAMA are in qualitative agreement with the 
experiment in all studied regions, i.e. in the large and low x limits as well as in the resonance region. 
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1 Introduction 



This work is devoted to modeling of the scattering amplitude for inelastic electron-proton scattering. 
The kinematics of inclusive ep scattering, applicable to both high energies, typical of HERA, and low 
energies as at JLab, is shown in Fig. ^ We introduce virtuality Q^, = —q^ = —{k — k'Y > 0, and 
'sj" ' Bjorken variable x = Q^/2p ■ q. These variables x, and Mandelstam variable s (of the 7*p system), 
I s — (p + qY, obey the relation: 

s = Q^{l-x)/x + m'^ , (1) 
where m is the proton mass. And Fig. |21 shows how inelastic 7*p scattering is related to the forward 



Ph. elastic (t=0) 7*p scattering, and then the latter is decomposed into a sum of the s— channel resonance 



exchanges. 
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Figure 1: Kinematics of inelastic electron-proton scattering. 

About thirty years ago Bloom and Gilman 5 observed that the prominent resonances in inelastic 
e~p scattering (see Fig. ^ do not disappear with increasing photon virtuality , but fall at roughly the 
same rate as background. Furthermore, the smooth scaling limit proved to be an accurate average over 
resonance bumps seen at lower and s, this is so called Bloom-Gilman or hadron-parton duality. Since 
the discovery, the hadron-parton duality was studied in a number of papers |S] and the new supporting 
data has come from the recent experiments [3 These studies were aimed mainly to answer the 
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questions: in which way a limited number of resonances can reproduce the smooth scahng behaviour? The 
main theoretical tools in these studies were finite energy sum rules and perturbative QCD calculations, 
whenever applicable. Our aim instead is the construction of an explicit dual model combining direct 
channel resonances, Regge behaviour typical for hadrons and scaling behaviour typical for the partonic 
picture. Some attempts in this direction have already been done in Refs. 0121210, which we will discuss 
in more details below. 
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Veneziano duality 

Figure 2: According to the Veneziano (or resonance-reggeon) duality a proper sum of either t-channel or 
s-channel resonance exchanges accounts for the whole amplitude. From 0. 



The possibility that a limited (small) number of resonances can build up the smooth Regge behaviour 
was demonstrated by means of finite energy sum rules [Sj. Later it was confused by the presence of 
an infinite number of narrow resonances in the Veneziano model [TO , which made its phenomenological 
application difficult, if not impossible. Similar to the case of the resonance-reggeon duality 9 , the 
hadron-parton duality was established 5 by means of the finite energy sum rules, but it was not realized 
explicitly like the Veneziano model (or its further modifications). 

First attempts to combine resonance (Regge) behaviour with Bjorken scaling were made ^2 El El at 
low energies (large x), with the emphasis on the right choice of the Q^-dependence, such as to satisfy the 
required behaviour of form factors, vector meson dominance (the validity (or failure) of the (generalized) 
vector meson dominance is still disputable) with the requirement of Bjorken scaling. Similar attempts 
in the high-energy (low x) region became popular recently stimulated by the HERA data. These are 
discussed in section O 

Recently in a series of papers 0(2112101 authors made attempts to build a generalized Q^-dependent 
dual amplitude A{s,t,Q^). This amplitude, a function of three variables, should have correct known 
limits, i.e. it should reduce to the on shell hadronic scattering amplitude on mass shell, and to the nuclear 
structure function (SF) when t = 0. In such a way we could complete a unified " two-dimensionally dual" 
picture of strong interaction 0|21I210 - see Fig. O 

In Ref. the authors tried to introduce Q^-dependence in Veneziano amplitude or more 
advanced Dual Amplitude with Mandelstam Analyticity (DAM A) ^11. The Q^-dependence can be intro- 
duced either through a Q^-dependent Regge trajectory 0, leading to a problem of physical interpretation 
of such an object, or through the g parameter of DAMA 00. This last way seems to be more realistic 
, but it is allowed only in the limited range of due to the DAMA model requirement g > 1 (see 
for details). 

In the papers the authors went in an opposite direction - they built a Regge-dual model with 
Q^-dependent form factors, inspired by the pole series expansion of DAMA, which fits the SF data in 
the resonance region^. The hope was to reconstruct later the Q^-dependent dual amplitude, which would 
lead to such an expansion. 

A consistent treatment of the problem requires the account for the spin dependence, which we ignore 
in this paper for the sake of simplicity. Our goal is rather to check qualitatively the proposed new way 
of constructing the "two-dimensionally dual" amplitude. 

^It is important that DAMA not only allows, but rather requires nonUnear complex Regge trajectories |14|. Then the 
trajectory with restricted real part lead to a limited number of resonances. 
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Figure 3: Veneziano, or resonance-reggeon duality 10 and Bloom-Gilman, or hadron-parton duality |S] 
in strong interactions. From 12]. 



2 Modified DAMA model 



The DAMA integral is a generalization of the integral representation of the B-function used in the 
Veneziano model 23| ^: 

where a' = a(l — z), a" = az, and 5 is a free parameter, g > I, and as{s) and at{t) stand for the Regge 
trajectories in the s— and channels. 

In this paper we propose a modified definition of DAMA (M-DAMA) with Q^-dependence ^Hl- It 
also can be considered as a next step in generalization of the Veneziano model. M-DAMA preserves all 
the attractive features of DAMA, such as pole decompositions in s and t, Regge asymptotics etc., yet it 
gains the Q^-dependent form factors, correct large and low x behaviour for i = etc. 

The proposed M-DAMA integral reads |15) : 

(^) . (3) 

where /3((5^) is a smooth dimensionless function of Q^, which will be specified later on from studying 
different regimes of the above integral. 

The on mass shell limit, = 0, leads to the shift of the s— and i— channel trajectories by a constant 
factor /3(0) (to be determined later), which can be simply absorbed by the trajectories and, thus, M- 
DAMA reduces to DAMA. In the general case of the virtual particle with mass M we have to replace 
by ((52 + M2) in the M-DAMA integral. 

^There are several integral representations of DAMA 1141 . here we shall use the most common one. 
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Now all the machinery developed for the DAM A model (see for example ^1]) can be applied to the 
M-DAMA integral. Below we shall report briefly only some of its properties, relevant for the further 
discussion. 



2.1 Singularities in M-DAMA 



The dual amplitude D(s, t, Q^) is defined by the integral Q in the domain TZe (as(s') < and 

TZe {at{t") + I3{Q'^')) < 0. For monotonically decreasing function TZe l3{Q'^) (or non-monotonic function 
with maximum at — 0) and for increasing or constant real parts of the trajectories these equations, 
appHed for < z < 1, mean TZe (ctsis) + /3(0)) < and TZe {at{t) + /3{0j) < . To enable us to study 
the properties of M-DAMA in the domains TZe (a^(s') + f3{Q'^")) > and TZe {at{t") + /3(Q^')) > 0, 
which are of the main interest, we have to make an analytical continuation of M-DAMA. This leads to 
the appearance of two moving poles 

as(s(l - z„)) + /3(Q^z„) = n and 
at{tz„^)+l3{Q^{l-z„,))=m, ti, m = 0, 1, 2... (4) 

The singularities of the dual amplitude are generated by pinches which occur in the collisions of the above 
mentioned moving and fixed singularities of the integrand z = 0, 1. 



1 . The collision of a moving pole z = z„ with the branch point z = results in a pole at s = s„ , where 
s„ is defined by 

a,(s„)-f /3(0) =n. (5) 

Please, notice the presence of an extra (in comparison to DAMA) term /3(0). It can be considered 
as a shift of the trajectory. If /3(0) is an integer number, then the modification is trivial. 

2. The collision of a moving pole z = z„ with the branch point z — 1 results in a pole at — Q^, 
defined by 

asiO)+(3{Ql)^n. (6) 

In this sense we can think about (3{Q^) as of a kind of trajectory, but we do not mean that it 
describes real physical particles. Also we will see later that with a proper choice of (3{Q^) we can 
avoid these unphysical poles, and P{Q'^) required by the low x behaviour of the nucleon SF is exactly 
of this type. 

3. Similarly, the collision of a moving pole z = z„j with the branch point z = 1 results in a pole at 
t = tm, defined by 

at{tm)+m)=m. (7) 

4. The collision of a moving pole z = Zm with the branch point z = results in a pole at ~ Q^m 
defined by 

am + ml^)^rn. (8) 

Note that if as(0) — at(0) the poles in will be degenerate. For further discussion we shall 
consider a non-degenerated case. 



2.2 Pole decompositions 



Similarly as for DAMA ^3], case 1 from the above results into pole decomposition of M-DAMA amplitude 
with the following expression for the pole term |15| : 
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where 



Ci{t, 



1 rf' 



1 - z 



-at{tz)-p{Q^(l^z))-l 



(10) 



z=0 



Formula @ shows that our D{s,t,Q'^) does not contain ancestors and that an (n + l)-fold pole emerge 
on the n-th level. The crossing-symmetric term can be obtained in a similar way by considering the case 
3 from the list above. 

The modifications with respect to DAMA are A) the shift of the trajectory as{s) by the constant 
factor of /3(0) (we can easily remove this shift including /?(0) into trajectory); B) the coefficients C; are 
now Q^-dependent and can be directly associated with the form factors. The presence of the multipoles, 
eq. ijnji, does not contradict the theoretical postulates. On the other hand, they can be removed without 
any harm to the dual model by means the so-called Van der Corput neutralizer^ . This procedure ^1] 
seems to work for M-DAMA equally well as for DAMA and will result in a " Veneziano-like" pole structure: 



C„(t, Q2) 



n - a,{s) - I3{0) 



(11) 



The Q^-pole terms can be obtained by considering cases 2 and 4 from section ITTl but, as we shall see 
later in section^ with our choice of P{Q^) we avoid poles. 



2.3 Asymptotic properties of M-DAMA 

Let us now discuss the asymptotic properties of M-DAMA. Using exactly the same method as in H^H it is 
possible to show that if the trajectory satisfies some restriction on its increase, then we obtain the Regge 
asymptotic behaviour |15| : 

i?(s,t,Q2)^s"*W+'3(0)/('3'), s->oo. (12) 

So, in the Regge limit M-DAMA has the same asymptotic behaviour as DAMA (except for the shift 
/3(0)). 

It is more interesting to study the new regime, which does not exist in DAMA - the limit — s- oo, 
with constant s, t. We assume that l3{Q'^) — > — oo for — > oo. Then |15| . 

where W « 8j\n{Q'^/QQ). For DIS, as we shall see below, if s and t are fixed and oo then 

u — —2Q^ — ^ — oo, as it follows from the kinematic relation s + 1 + u — 2m? — 2Q^. So, we need also 
to study the D{u,t,Q^) term in this limit. If |q;„(— 2Q^)| is growing slower than |/3((5^)| or terminates 
when oo, then the previous result (eq. (|13|l . s to be changed to it = —2Q^) is still valid. 



3 Nucleon structure function 



The total cross section of 7*p scattering is related to the SF by 



^In brief, the procedure 1141 is to multiply the integrand of J^l t>y a function 4>{z), which has the following propcrti( 

0(0) =0, 0(1) = 1, 0"(1) = O, n = 1,2,3,... 



The function 4){z) = 1 — exp — | , for example, satisfies the above conditions. 
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where a is the fine structure constant. In cq. ifT^ we neglected R(x, Q^) = cfl{x, Q'^)/(Tt{x, Q^), which 
is a reasonable approximation. 

The total cross section is related to the imaginary part of the scattering amplitude 

af^ix, Q') = Im Q^), i = 0, Q^) . (15) 



where Pcm is the center of mass momentum of the reaction, Pcm = 2(i-x) y for DIS. Thus, 



P2 {x, Q') = ,w/, . ,;^,^v. ^(*(^' Q^),t = Q,Q^). (16) 



we have 

a (s - m^) (1 + 4m2a;2/(32)3/2 
The minimal model for the scattering amplitude is a sum |17| 

A{s, 0, Q2) = c(s - u){D{s, 0, Q2) _ 0, Q")), (17) 

providing the correct signature at high-energy limit, where c is a normalization coefficient. As it was said 
at the beginning, we disregard the symmetry properties of the problem (spin and isospin) , concentrating 
on its dynamics. 



In the low x limit: a; — > 0, i = 0, = const, s = Q"^ jx — > cxd, u ~ ~s we obtain from eqs. H12I17(I : 

Im Ais,0,Q^)\s-,oo ~ s^'W+z^W+i/W') . (18) 



Our philosophy in this section is the following: we specify a particular choice of /3(Q^) in the low x 
limit and then we use M-DAMA integral lO to calculate the dual amplitude, and correspondingly SF, 
in all kinematical domains. We will see that the resulting SF has qualitatively correct behaviour in all 
regions. Even more - our choice of (3{Q'^) will automatically remove poles. 

According to the two-component duality picture |18| , both the scattering amplitude A and the struc- 
ture function f 2 sue the sums of the diffractive and non-diffractive terms. At high energies both terms are 
of the Regge type. For scattering only the positive-signature exchanges are allowed. The dominant 
ones are the Pomeron and / Reggeon, respectively. The relevant scattering amplitude is as follows: 

— ' (19) 
/ 

where at and Rk are Regge trajectories and residues and k stands either for the Pomeron or for the 
Reggeon. The residue is chosen to satisfy approximate Bjorken scaling for the SF |19l I2U| . From eqs. 
H16I19|) SF is given as: 

(Q \ Qfc (0) —1 
— ■ (20) 

Bjorken variable x — Q'^ /s for s ^ 00 and thus, Regge asymptotics and scaling behaviour require that 

RkiQ') - (QT'"''^"^ ■ (21) 

Actually, it could be more involved if we require the correct ^ limit to be respected and the 
observed scaling violation (the "HERA effect") to be included. Various models to cope with the above 
requirements have been suggested [SlEIEni- At HERA, especially at large Q^, scaling is so badly 
violated that it may not be explicit anymore. 

In the phenomenological models which are used nowadays to fit F2 data |19[ I2(JI [71 |H1 1^ (also 0] 
were discussed in introduction) the Q^-dependence is introduced " by hands" , via residue in the form (|21|l , 
parameters of which are then fitted to the data. Now we have a model which contains Q^-dependence 
from the very beginning and automatically gives a correct behaviour of the residues. 
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Data show that the Pomeron exchange leads to a rising structure function at large s (low x). To 
provide for this we have two options: either to assume supercritical Pomeron with ap{0) > 1 or to 
assume a critical (ap(0) ~ 1) dipole (or higher multipole) Pomeron 1211 E^ . The latter leads to the 
logarithmic behaviour of the SF: 

J^2,p(x, ^ Q2^p(g2) ln(-^) , (22) 

which proves to be equally efficient |16[ ■ 

Let us now come back to M-DAMA results. Using eqs. (|16I18II we obtain: 

F2 ~ s"*(0)+/3(0)q2^^(Q=) _ (23) 

Choosing 

PiO) = -1 (24) 

we restore the asymptotics (l2Uf) and this allows us to use trajectories in their commonly used form. Now 
we have to find such a PiQ"^), which can provide for Bjorken scaling. If we choose (3{Q'^) in the form 

/3(g2) =d-7ln(gVQ2)^ (25) 

with 

7 = («t(0) + /3(0) + l)/lng = atiO)/lng, (26) 
where d, Qq are some parameters, we get the exact Bjorken scaling. 

Actually, the expression (|25|l might cause problems in the ^ limit. To avoid this, it is better to 
use a modified expressions 



This choice leads to 

F.(x,<5V..-.<».(^) , (28, 

where slowly varying factor (^ Q2^ga 1 is typical for the Bjorken scaling violation (for example j2()|). 



Now let us turn to the large x limit. In this regime a; ^ 1, s is fixed, = ~^ °° ^^'^ 

correspondingly u = —2Q^. Using eqs. H13I16I17|I we obtain: 

^2 ~ (1 - x)2gV''(^'/'^ (g"4V2) „ ^a„(-Q=)^ . (29) 

For Q"^ ^ oo factors — ^""(-"3^)^ and W are slowly varying functions of under our assump- 

tion about q;m(— g^). Thus, we end up with qualitatively correct behaviour 

F^^i^j _ (i_^)2o,(0)ln2<,/ln5^ (gg^ 

Let us now study F2 given by M-DAMA in the resonance region. The existence of resonances in SF 
at large x is not surprising by itself: as it follows from H15|l and (|16|l they are the same as in 7*73 total 
cross section, but in a different coordinate system. 

For M-DAMA the resonances in s-channel are defined by the condition ©. For simplicity let us 
assume that we performed the Van der Corput neutralization and, thus, the pole terms appear in the 
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form In the vicinity of the resonance s — sj^^g only the resonance term Dj^es(s, 0, Q^) is important 
in the scattering ampHtude and correspondingly in the SF. 



Using (3{Q'^) in the form l|27|l . which gives Bjorken scaling at large s, we obtain from eq. H10|l : 



, "t(0) 
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Q2 + Q2 Ij,^^ ^ Q2 



^0 



(31) 



The term ( q^+q-^ ) gives the typical Q -dependence for the form factor (the rest is a slowly varying 
function of Q^). 

If we calculate higher orders of C„ for subleading resonances, we will see that the Q^-dependence is 

/ q2 \at(0) 

still defined by the same factor ( Q2^Qi I • Here comes the important difference from the Regge-dual 

model 131 13] motivated by introducing Q^-dependence through the parameter g. As we see from eq. (|ll|l . 
g enters with different powers for different resonances on one trajectory - the powers are increasing with 

the step 2. Thus, if g ~ ^ q^+q'^ ) ' ^^^^ form factor for the first resonance is [n — Q) ^ (q^^Q^) ' 

and for the second one {n = 2) it is ^ \ W+q^) discussed in 0) the present accuracy of 

the data does not allow to discriminate between the constant powers of form factor (for example Refs. 
123 El El El) and this work) and increasing ones. 



4 How to avoid poles? 

General study of the M-DAMA integral allows the existence of poles (see cases 2, 4 in section |OJ 
which would be unphysical. The appearance and properties of these singularities depend on the particular 
choice of the function /3((5^), and for our choice, given by eq. 1)27(1 . the poles can be avoided. 

We have chosen /3(Q^) to be a decreasing function, then, according to conditions (|6I8|I . there are no 
poles in M-DAMA in the physical domain > 0, if 

Tie ^(0) < -a^(0) , Tie ^(0) < -at(0) . (32) 

We have already fixed /3(0) = —1, eq. (|24|l . and, thus, we see that indeed we do not have poles, except 
for the case of supercritical Pomeron with the intercept ap(0) > 1. Such a supercritical Pomeron would 

-' pole 



generate one unphysical pole at = Qvoie defined by equation 



i_^in(^^)+c.p(o) = o ^ QLie-Qlig^^-i)- (33) 



ln.g V Qo 

Therefore we can conclude that M-DAMA does not allow a supercritical trajectory - what is good property 
from the theoretical point of view, since such a trajectory violates the Froissart-Martin limit j25| . 

As it was discussed above there are other phenomenological models which use dipole Pomeron with 
the intercept ap{0) = 1 and also fit the data (see for example JHD- This is a very interesting case - 
(at{Q) = 1) - for the proposed model. At the first glance it seems that we should anyway have a pole at 

= 0. It should result from the collision of the moving pole z = zq with the branch point z = 0, where 
at(0) + /3{Q^{1 — zq)) = in our case. Then, checking the conditions for such a collision: 

atiO) -ta'M^o + - fi'{Q^)Q^zo = ^ zo = ""'^'^^ ~ ^^^^^ 



we see that for t = and for /3(Q^) given by eq. (|27|l the collision is simply impossible, because zo(Q^) 
does not tend to for 0. Thus, for the Pomeron with ap{0) = 1 M-DAMA does not contain any 

unphysical singularity. 
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On the other hand, a Pomeron trajectory with ap(0) = 1 does not produce rising SF as required 
by the experiment. So, we need a harder singularity and the simplest one is a dipole Pomeron. A 
dipole Pomeron produces poles of the second power - Ddipoie{s,tm) oc (^m-ap{t)+iy^ ' example 
ref. \21\ and references therein. Formally such a dipole Pomeron can be written as — 



dap (m — ap + ' 

and generalizing this - Ddipoie{s, t) — g^D{s, t) , where D{s, t) can be given for example by DAMA or 
M-DAMA. Applying this expression to the asymptotic formula of M-DAMA, eq. H12|l . we obtain a term 
gPiQ ) gat(t)+p{o) i^g^ which then leads to a logarithmically rising SF (for ap(0) + /3(0) = 0) - the one 
given by eq. (|^ . 

For (3{Q'^) in the form H27() M-DAMA will generate an infinite number of the poles concentrated 
near the "ionization point" = —Qq. Although these are in the unphysical region of negative Q^, such 
a feature of the model 

A) makes us think about (3{Q'^) as about a kind of trajectory, what is not the case, as it was stressed 
above, and 

B) might create a problem for a general theoretical treatment, for example for making analytical con- 
tinuation in . To avoid this we can redefine l3(Cp) in the nonphysical region, for example in the 
following way: 

r _i_££lMinf«li^) , forQ2>o, 

This function has a maximum at = 0, /3(0) — —1. M-DAMA with /3((5^) given by eq. (|34|l preserves 
all its good properties, discussed above, and docs not contain any singularity in (except for the 
supercritical Pomeron case, which we do not allow). 



5 Conclusions 

A new model for the Q^-dependent dual amplitude with Mandelstam analyticity is proposed. The 
M-DAMA preserves all the attractive properties of DAMA, such as its pole structure and Regge asymp- 
totics, but it also leads to generalized dual amplitude A{s,t,Q'^) and in this way realizes a unified 
"two-dimensionally dual" picture of strong interaction El 13 (see Fig. This amphtude, when 

t = 0, can be related to the nuclear SF, and in this way we fix the function (i{Q^), which introduces the 
Q^-dependence in M-DAMA, eq. Our analyzes shows that for both large and low x limits as well as 
for the resonance region the results of M-DAMA are in qualitative agreement with the experiment. 

In the proposed formulation a Q^-dependence is introduced into DAMA through the additional func- 
tion /3((5^). Although in the integrand this function stands next to Regge trajectories, this, as it was 
stressed already, does not mean that it also corresponds to some physical particles. There is no qualita- 
tive difference between two ways of introducing Q^-dependence into DAMA: through the Q^-dependent 
parameter g, i.e. function g{Q^) EI or through the function (3{Q'^). On the other hand the second 
way, i.e. M-DAMA, is applicable for all range of and it results into physically correct behaviour in all 
tested limits. 
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